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Abstract 

We examine the quantum motion of two particles interacting through a con- 
tact force which are confined in a rectangular domain in two and three dimen- 
sions. When there is a difference in the mass scale of two particles, adiabatic 
separation of the fast and slow variables can be performed. Appearance of 
the Berry phase and magnetic flux is pointed out. The system is reduced 
to a one-particle Aharonov-Bohm billiard in two-dimensional case. In three 
dimension, the problem effectively becomes the motion of a particle in the 

presence of closed flux string in a box billiard. 
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The quantum billiard problem has been instrumental in revealing non-trivial features 
of seemingly simple systems with minimal degrees of freedom. As a consequence, unlike 
real-world billiards, quantum physicists' version of the billiard problem has only one particle 
moving within a walled boundary. The dynamics is controlled by the shape of the boundary, 
and by optional obstacles placed inside the billiard |TJ-|5|]. The formulation of zero-size 
obstacle offers a possibility to go beyond this limitation, since it can be thought of as 
the second particle with infinitely heavy mass. 

In this Letter, we consider the quantum mechanics of two and three dimensional systems 
consisting of two particles, one heavier than the other, interacting through a contact force, 
and confined in a walled rectangular boundary. Adopting the adiabatic separation of fast 
motion of the light particle and slow motion of the heavy one, we can formulate the problem 
in two steps. Namely, we first solve the motion of the light particle with the location of the 
heavy particle fixed, then calculate the effective Hamiltonian for the heavy particle using 
the eigenstates of the light particle. The first step yields a one-particle billiard problem 
in the presence of a pointlike scatterer, which has been the focus of several recent studies 
||7| PH[. The identification of the Berry phase [L]]] in this one-particle billiard |9| then has 
an interesting ramification in the second step: The effective Hamiltonian for the heavy 
particle acquires a non-trivial vector potential which results in a field formally identical to 
the magnetism. We show that, in two dimension, this "magnetic" field appears as infinitely 
thin flux lines penetrating the billiard. Through a numerical example, we demonstrate that 
the flux lines have discernible effects both on the spectra and wave functions of the system. 
We also show that, in three dimension, the flux lines form closed strings which reside inside 
the billiard. Several examples of the flux string are displayed for a box billiard. 

Let us suppose that two particles m and M each located at f and R interacting through 
a delta force are confined in a d- dimensional rectangular domain by a hard wall V wa u. We 
assume that M is substantially larger than m so that the motion of M can be treated 
adiabatically. The Hamiltonian describing the system is given by 
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where h(R) is a sub-Hamiltonian for the light particle which is given by 



-Q 

V 



(2) 



h(R) = ^- + v6 d (R - f) + V wall {r) 
2m 



For a fixed value of R, this sub-system is known as the billiard with a pointlike scatterer 
As it stands, the resolvent of this sub-Hamiltonian is not well-defined for spatial dimension 
d > 2. However, in case of d < 3, the system can be made meaningful by redefining a 
coupling constant v in terms of renormalized coupling constant v. Assuming that the h(R) 
is handled with such proper renormalization procedure || , we write the eigenvalue equation 
of h(R) as 



h{R) a(R)) = u a {R) a{R)) . (3) 

There are two types of solutions to this equation ||. Let us denote the eigenvalue and 
eigenstate with v — (empty billiard) as e n and |n). The first, trivial type of solution exists 
when the delta potential is located on the node-line of unperturbed state; 



u a = e n if f n {R) = 0. 



(4) 



where / n (^) = represents the node-line of the wavefunction (r| n). The second, more 
generic type of solution of eq. (0) is obtained from the equation 



G(R;u) - - = 
v 



(5) 



where 



G(r] uj) = (r| n){n\ r) 



UJ - £ n 



+ 



e 2 + 1 



(6) 



is the inverse of transition matrix of the light particle m. If we assume that m stays at 
the adiabatic state 



a(R)) during the motion of M, we obtain an effective Hamiltonian 



governing the dynamics of M as 



2M 

where the scalar potential U a is given by 



Hf = (a(R) H a{R)) (7) 

\t(p- MR)) 2 + u a (R) + V wall (R) 



U a (R)=u a (R), (8) 



and the gauge potential (Mead-Berry connection) A a is given by [[Tll ,|l2 



A a {R)=i(a{R) V R a(R)). (9) 



It is important that we choose normalized eigenstate a(R)j as a single- valued function 
in the parameter space of R. In some cases, A a can be made to be zero everywhere by 
redefining the phase of the a(R)^, but in general, this is not possible. When R makes an 
adiabatic cyclic motion, it acquires a geometric phase given by the circular line integral of 
A a . The appearance of the vector potential is a consequence of the fact that this phase 
is non-zero. It is worth stressing that in the system considered here, the parameter space 
of slow variable is nothing but the coordinate space of heavy particle R. Thus the vector 
potential A a resides in the coordinate space R. The eigenvalue equation 

Hf\$ Qp ) = E ap \$ ap ) (10) 

gives the energy states of the whole system E ap along with the eigenstate of the heavy 
particle They are indexed by two integers a and p each specifying the state for light 

and heavy particles. 

We first concentrate on the case of two dimension. Since the Hamiltonian is real sym- 
metric, the geometric phase is tt when the path of the cyclic motion of R surrounds the 
diabolical point R* which is the point of degeneracy of the sub-Hamiltonian h(R); 

namely u> a (R*) = u a ±i(R*). Otherwise, the phase is zero. Since the quantity G(f; uj) is a 
monotonous function of uj except at it's poles M, the degeneracy can occur only between 
the solutions of eqs. @ and (5). Therefore, location of R* is determined by 
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f n (R*) = and g n (R*) = (11) 

where g n (R) is defined by 

g n (R) = G(R;£ n )-~. (12) 

In general, this equation has more than one solution for each degeneracy e n = u a . We 
write them as R* a j (j = 1,2, ■ ■ •). Because of the sign reversions around R^j, single- valued 
representation of the eigenfunction a(i?)y necessarily becomes complex even for the real 
symmetric Hamiltonian. There is a freedom to choose the gauge. We take the "maximally 
symmetric" choice around each i2* ■; 

z„ - „ (n\R 



a(R)) = exp (- £ ip aj )N a (R) £ |n) ) 1 (13) 

where ip a j is an angle variable of R in the polar coordinate with its origin at i2* • , and N a is 
a real normalization factor. From eqs. <M) and (p!3|), we obtain 

Ax{R) = e - r\ ( i4 ) 



. 2 



where e Vaj is a unit vector of the direction ip a j. We choose its direction such that (R—R^j) x 
e ¥ , Qj . points to the positive z-axis in usual right-hand convention. The gauge curvature (or 
the "magnetic" field) B a is given by 

B a (R) = V x A a (R) = J2(-nz)*8 2 (R ~ K 3 ) (15) 

i 

where n 2 is the unit vector of the direction of z-axis. The problem is now turned into 
the motion of a particle M in the presence of infinitely thin lines of magnetic flux piercing 
through the billiard. This type of system is known as Aharonov-Bohm billiard |H,I5 



Each flux line has the strength 7T, or 1/2 in the unit of 2tt which is usually adopted in the 
literature. This is known to give the anti-unitary symmetry to the system guaranteeing 
the correct global quantum properties such as level statistics. It is easy to see that the 
eigenvalues of the effective Hamiltonian, eq. (0) are unchanged under n z — > —n z , while the 



phase of eigenf unctions changes. This means that the physics of a particle M is invariant 
under the reversion of the vector potential. In fact, for infinitely thin flux, the direction of B 
can be changed arbitrarily as long as it is off the billiard plane. This corresponds to the gauge 
transformation of the vector potential. Explicit calculation of eq.(|7]) with eq.flHj) shows that 
the vector potential acts as centrifugal barrier with square inverse radial dependence around 
each diabolical points. That will effectively cause the wavefunctions to avoid the flux lines. 

To illustrate our arguments, we show a numerical example. The billiard boundary is 
chosen to be a rectangle of size [L x 1/L] where L = 1.1557. The mass of the light particle 
is set to m = 2tt, and the heavy particle M = 5m. The strength of the coupling is chosen 
to be v = 10. In Fig. 1, the energy eigenvalue of the system is shown. In each of the energy 
levels, left side is the full calculation, and the right side is the calculation neglecting the 
effects of flux lines. The density profile of the eigenstates of the heavy particle (R |$ ap ) 
for selective a and p are depicted in Fig. 2, along with the location of the magnetic flux 
which is shown as black squares. The repulsive effect of the magnetic flux line is clearly 
visible. 

We now turn to the three dimensional case. We consider a generic case of three- 
dimensional box-shaped domain in which unperturbed single-particle spectra have no degen- 
eracies. The diabolical locations are determined by the same eq. (|ll|) as in two dimensional 
case, but with R* now representing a three dimensional coordinate vector. The condition, 
eq. (|TT|) specifies a one- dimensional line in the space of R. This is to be expected also from 
another perspective. Since the degeneracy is co-dimension two for real symmetric Hamil- 
tonian, R* cannot be the isolated points in the three-dimensional space of R. Rather, it 
becomes a line of degeneracies. Clearly, the line cannot have the edge inside the billiard 
domain. Also, the line does not go through the domain wall, since this would mean that the 
unperturbed system has a degeneracy, which we have excluded in our assumption. The only 
possible configuration for the line of degeneracy, therefore, is the closed string. When one 
makes adiabatic circular motion of R, the state of the sub-system a(R)"j obtains a phase 
7r if the circle is such that degeneracy strings pierce through its area odd number of times. 



Otherwise, the geometric phase is zero. With the Stokes theorem, it is easy to show that 
the gauge potential A a whose circular line integral has such properties gives the flux B a 
which is non-zero only along the degeneracy string. Thus the string of degeneracy might be 
called the flux string. The strength of the flux B a along the string is 1/2 in the unit of 2n. 
Formally, it is given by 

B a {R) = Q / rf J R 5 5(/ n ( J R 5 ))5(^( J R 5 ))n( J R 5 )27r5 3 ( J R - R s ) (16) 

where n(R) is the direction vector of the flux string. As in dimension two, one can reverse 
the direction of n without changing the physics. 

As in two dimension, the flux string in three dimension will have observable effects on the 
motion of the slow particle, and thus on the energy spectra of the system. Such calculation 
tends to become tedious for three dimension. Instead, in Fig. 3, we simply display the flux 
string for the system of two particles in a box of size [L x ,L y , l/(L x L y )} where L x = 1.1557 
and L y = 1.0310. The mass and the coupling constant are chosen to be m = 2n and v = 10 
as before. It is evident that the flux strings possess the pictorial nature which is quite 
striking. We think that the impact is enhanced by the fact that these objects emerge from 
nothing more than two interacting particles in a box. If the observation of wave functions in 
a box becomes possible, these flux strings will be visible as the empty void of the probability 
distribution. 

The magnetic flux string similar to the one found here has been predicted and discussed 
in the context of field theories. In elementary particle physics, it is known as U(l) vortex 
T6| , and in cosmology as cosmic string ||17|| . It is quite conceivable that the actual exper- 



imental observation of magnetic flux string may be first realized at mesoscopic scale in a 
system similar to the one considered here. Further, the present study with a very simple set- 
ting suggests that, despite its exotic appearance, the flux string is ubiquitous in many-body 
quantum systems. We conclude that the quantum billiards acquire considerable richness 
with the introduction of more than one moving particle. 
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FIGURES 

FIG. 1. Energy level scheme of the eigenvalue equation eq. (10) with the effective Hamiltonian, 
eq. (7). The size of the rectangular boundary is set to [1.15572,0.86526]. Other parameters are 
m = 2tt, M = IOtt and v = 10. 

FIG. 2. The density profiles of the probability distribution of the heavy particle of the several 
selective eigenstates of eq. (10) for the rectangular billiard of L = 1.15572. the numbers in the 
figure refers to the quantum numbers (a,p). 

FIG. 3. Flux strings found in the two particle billiard problem in a three dimensional box of 
the size [1.1557, 1.0310,0.8392]. (a) n = 2, (b) n = 3, (c) n = 4, (d) n = 5, (e) n = 6, (f) n = 7, 
in eq. (11), respectively. 
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Fig. 1 
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